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ABSTRACT 

We present a method for obtaining accurate black hole (BH) mass estimates 
from the Mg II emission line in active galactic nuclei (AGNs). Employing the 
large database of AGN measurements from the Sloan Digital Sky Survey (SDSS) 
presented by Shen et al., we find that AGNs in the redshift range 0.3-0.9, for 
which a given object can have both H/3 and Mg II line widths measured, display 
a modest but correctable discrepancy in Mg Il-based masses that correlates with 
the Eddington ratio. We use the SDSS database to estimate the probability 
distribution of the true (i.e., H/3-based) mass given a measured Mg II line width. 
These probability distributions are then applied to the SDSS measurements from 
Shen et al. across the entire Mg Il-accessible redshift range (0.3-2.2). We find 
that accounting for this residual correlation generally increases the dispersion of 
Eddington ratios by a small factor (~0.09 dex for the redshift and luminosity bins 
we consider). We continue to find that the intrinsic distribution of Eddington 
ratios for luminous AGNs is extremely narrow, 0.3-0.4 dex, as demonstrated by 
Kollmeier et al. Using the method we describe, Mg II emission lines can be used 
with confidence to obtain BH mass estimates. 
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1. INTRODUCTION 

The Mbh-c" correlation between super- 
massive black hole (BH) mass and the ve- 
locity dispersion of the surrounding stel- 
lar system indicates a significant connec- 
tion between galaxy and BH assembly 
( Ferrarese fc Merritt 2000l : Gebhardt et al. 
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2000). Thus, it is essential to establish the 
most basic BH parameters: their intrinsic 
distribution of masses and growth rates. 

It is extraordinarily difficult to measure 
BH masses directly, even in the nearby uni- 
verse, because of the small spatial scales 
that must be resolved to probe the grav- 
itational influence of the BH. In active 
galactic nuclei (AGNs ), the technique of re- 
verberation mapping (IBlandford &: McKee 
19821 : iPetersoru Il993l ) employs high reso- 
lution in the time domain to probe gas 
dynamics on spatial scales close to the 
BH. However, the long-duration spectro- 
scopic monitoring campaigns required for 
reverberation studies currently preclude 
the method from being applied to large 
numbers of objects. Therefore, one must 
rely on even more indirect techniques of 
BH mass estimation to build up statisti- 
cally signiflcant samples. 

The "virial" method has been empiri- 
cally calibrated fro m reverberation map- 



ping experiment s (IWandel et al.l 11999 



Vestergaardll2002l : iMcLure fc Jarvid l2002t ) 



and allows a BH mass estimate from a 
one-time measurement of the width of a 
broad emission line and the AGN lumi- 
nosity (see § 2). To facilitate applica- 
tion of the virial technique to large op- 
tical AGN surveys, versions have been de- 
veloped for H/3 at low redshift, for the 
Mg II doublet near 2800A at intermedi- 
ate redshift, and for the C IV doublet 



prescriptions, see 


Vestereaard & Peterson 


2006; 


McGill et al. 


2008). 



In this paper, we present evidence of a 
systematic discrepancy in Mg Il-based BH 
mass estimates function of Eddington 
ratio (the ratio of the bolometric luminos- 



ity, Lboi, to the luminosity required for ra- 
diation pressure to balance the gravity of 
the BH), as well method to correct 

this trend. 

2. Method of Analysis 

The standard equation for estimating 
BH masses in AGNs from single-epoch 
spectroscopy is: 

log MBH = a + b log + 2 log y, (1) 

where V is the width of the broad emission 
line, Lc is the continuum luminosity near 
the line, and a and b are constants (which 
vary from line to line). The dependence on 
Lc arises because the distance of the broad 
emission line gas from the BH has been ob- 
served to correlate tightly with the AGN 
l uminosity oyer 4 o rders of magnitude in 



(IBentz et al.l 120071 ). Thus, the mass equa- 
tion reverts to the simple, virial combina- 
tion of radius and velocity. The calibra- 
tion of these relations rests on the bedrock 
of reverberation mapping measurements 
of local AGNs. By far, the best rever- 
beration mapping d ataset exists for H/3 



(jPeterson et al.ll2004l ). The a and b coeffi- 



cients for both C IV and Mg II principally 
rely on empirical correlations with B.f3- 
reverberation masses. However, whereas 
C IV reverberation studies of a handful 
of objects are consistent with e xpectations 



from B.p (jPeterson et al.ll2005l ). a clear in- 



dication of reverberation has not yet been 
found for Mg II. Although the line has been 



IClavel et al. 


1991; 


Dietrich & Kollatschnvl 


11995 




Metzroth et al. 




2006: 


WoQ .2008lj. 



only weak Mg II reverberat i on sig nals have 
been seen ( iReichert et al.l |l99J). Thus, 
the mass relation for Mg II relies on the 



2 



correlation of single-epoch estimates with 
B.p masses and the argument that Mg II 
and H/3 have similar ion ization potentials 
(jMcLure fc DunlopI 120041 ) . 

With the wide wavelength coverage of 
Sloan Digital Sky Survey (SDSS) spectra, 
certain redshift windows allow for two of 
these three emission lines to be measured 
simultaneously: B.f5 and Mg II are both 
accessible for z ~0.3-0.9, while Mg II and 
C IV can both be measured for z ~1.7- 
2.2. By comparing the BH mass estimates 
from the two lines, we are thus able to 
study systematic trends. Such a compari- 
son has been done previously for the mean 
relation, but the large sample size of the 
SDSS permits an analysis of higher-order 
correlations, which prove to be quite im- 
portant. 

Shen et al. (2008) have provided line 
width (FWHM) measurements and BH 
mass estimates — using the relations of 
McLure &: DunlopI (2004) and Vestergaard 
& Peterson (2006)— for roughly 60,000 
AGNs from the SDSS, including ~8,000 
with both H/3 and Mg II, and ~15,000 
with both Mg II and C IV. They provide 
a detailed analysis of the relationship be- 
tween Mg II and C IV, but simply describe 
the ratio of H/? to Mg II FWHM as follow- 
ing a log-normal distribution with a mean 
of 0.0062 dex and a dispersion of 0.11 dex. 
Shen et al. noted that the relation of the 
two FWHMs deviates slightly from a per- 
fect correlation, but did not explore the 
issue further. 

Under the premise that H/?, being the 
most extensively reverberation-mapped 
emission line, provides the best indicator 
of the BH mass, we examine AGNs having 
both B.(3 and Mg II mass measurements. 



Due to uncertainties in the line width mea- 
surements, we exclude the small number of 
AGNs that were flagged by Shen et al. as 
broad absorption line objects. However, 
the inclusion of these objects has negli- 
gible effect on our results. In Figure 1, 
we plot the difference in the (log of the) 
BH mass as a function of Eddington ratiqj. 
The strong correlation implies that if Mg II 
is calibrated simply from the mean of the 
H/?-Mg II relation, then it will underesti- 
mate the BH mass at low Eddington ratio 
and overestimate the mass at high Edding- 
ton raticl§. This would lead one to infer a 
narrower distribution of Eddington ratios 
than actually exists. One way to quantify 
this effect is to measure the slope of the 
correlation 



/3 = tan 



1 2Ci2 

- arctan 

2 C11-C22J 



(2) 



where Cij is the covariance matrix of the 
distribution shown in Figure [T] A slope 
(3 = would imply that Mg II provided 
completely independent information on the 
BH mass, while /3 = 1 would imply that no 
information is conveyed by the Mg II mea- 
surement. The actual value is (3 = 0.76, 
which means that Mg II is indicative of the 
BH mass, but must be treated with care. 

Simply modifying the luminosity- dependence 
of the Mg II mass formula (i.e., changing 
b in eq. [1]) cannot remove the observed 
trend. Therefore, to make a statistical 
correction to the Mg Il-based masses, we 
adopt the following method. For the SDSS 



^The Eddington ratio is computed with the H/3- 
based BH mass. 

^We note that the same trend from Figured] is seen 
when replacing FWHM with the inter-percentile 
value measurements of Fine et al. (2008). 
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objects with both R/3 and Mg II FWHMs, 
we look at AGNs with Mg II hnes in a 
given 0.1 dex bin of FWHM and tabulate 
the distribution of H/? FWHMs for those 
objects. We take those H/? distributions 
(normalized appropriately) as the proba- 
bility distributions for the true FWHM un- 
derlying the observed Mg II value (Fig. 2; 
Table 1). 

For any object with an accessible Mg II 
line, each bin in true FWHM is combined 
with the observed continuum luminosity 
to calculate a BH mass (via eq. [1]), and 
the mass is then used with the object's 
bolometric luminosity to derive an Edding- 
ton ratio. The probability for each bin of 
FWHM is added to the total number of 
objects contributing to the corresponding 
Eddington ratio bin (which are 0.2 dex- 
wide, allowing a one-to-one match between 
FWHM and Eddington ratio bins). Thus, 
each Mg II FWHM becomes a weighted 
distribution of Eddington ratio^, while 
H/9- and C IV-based Eddington ratios con- 
tribute directly at their observed values. 

We now apply this technique to the 
Shen et al. measurements of the uni- 
formly selected subsample of SDSS AGNs 
(Richards et al. 2006). 

3. Results & Discussion 

We construct distributions of Edding- 
ton ratios for SDSS in ranges of (Lboi, z) 
in Figure 3, showing the results both with 
(solid) and without (dotted) the Mg II sub- 



stitution presented above. The statistics 
of the distributions are given in Table 2. 
After correcting the Mg II measurements 
in the uniformly selected SDSS sample, we 
find that the average width of the Edding- 
ton ratio distribution increased by 0.09 dex 
for objects in the redshift range 0.3-2.2. 
We therefore find that the distribution of 
Eddington ratios remains very narrow at 



0.4 d ex, as was found by iKoUmeier et al. 



3 If the observed Mg II FWHM falls in a bin in 
which there were no H/3+Mg II measurements, it 
is given a probability of 1 within the bin corre- 
sponding to the Mg II FWHM. In the sample we 
consider, this applies to a single AGN. 



(120061) for the AGES-I survey. 

To test our procedure, we apply it to 
multiple subsamples for which we have 
both H/3 and Mg II data. Figure |4] shows 
the distribution of Eddington ratios for 
each subsample calculated in two ways, 
first using the true H/3 mass (solid) and 
second using our procedure applied to the 
Mg II mass (dashed). The similarity of 
these distributions (and the differences 
from the raw Mg Il-based values, shown 
as the dotted histograms) demonstrates 
that our procedure works well, recovering 
the true Eddington ratio distribution from 
the Mg II derived masses. 

The general trend we observe could be 
explained physically if the location where 
lines are formed in the broad-line region 
depends on accretion rate. In this case, 
the radius-luminosity relation would also 
depend on Eddington ratio and that would 
introduce an additional term in the virial 
mass relation. It would then be possible, 
in principle, to remove this dependence en- 
tirely analytically by fitting the observed 
correlation. We investigate this further in 
an upcoming work. 

The virial method for BH mass estima- 
tion has opened a new window in the study 
of supermassive BH demographics. While 
each indicator has systematics that must 
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be addressed (i.e., asymmetric line profiles, 
contamination from disk winds and metal 
lines, etc.), it is important to understand 
and attempt to correct for these systemat- 
ics so that this technique can be applied 
with confidence. In this contribution, we 
have identified a limitation in estimating 
BH masses from the Mg II line and pre- 
sented a way to remove it by exploiting 
the overlap of Mg II and B.f5 measurements 
presented by Shen et al. (2008). This 
method puts Mg II masses more securely 
on the same scale as H/3, which should be 
the most reliable in these studies. As the 
bias in the Mg II masses is also seen in our 
analysis of the 2dF Quasar Redshift Survey 
(2QZ) and the AGN and Galaxy Evolution 
Survey (AGES), in forthcoming work, we 
will apply our correction technique to those 
datasets (including the expanded sample of 
AGES-II). 

Understanding the transition between 
high-luminosity, high-redshift AGNs hav- 
ing a narrow Eddington ratio distribu- 
tion (~0.4 dex) and low-luminosity, low- 
redshift AGNs with a broad Eddington 
ratio distribution (>1 dex; Ho 2002; Woo 
& Urry 2002) can provide important con- 
straints on the physics of BH accretion. To 
determine these distributions, it is critical 
to continue to improve BH mass estimates 
that can be used through the bulk of the 
cosmic AGN activity. 
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Table 1 
Mg II Replacement Matrix 
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Note. — The values listed are the probabilities of the H/3 FWHM falling within a particular bin, given an input Mgll FWHM. All velocities are in log km s , and all 
bins are ±0.05 dex. 



Table 2 

Statistics of Eddington Ratio Distributions for SDSS 



Raw Mg II Replaced 
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3.27 



Note. — For each bin in redshift ^bin and luminosity ibin, the table lists: N, the number 
of AGNs; fx, the mean logarithm of the Eddington ratio; a, the dispersion in log(Lboi/^Edd); 
Sk, the skewness; and Ak, the kurtosis. 
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•..^g^?r^ajs.-:',.v:. 
■■■■J'iV.-.-- 



log (kol/LjIM) = log k.l 



-1 
[38.1 + log M(H/S)] 



Fig. 1. — Difference in H/? and Mg II esti- 
mates of log Mbh as a function of tlie log of 
the Eddington ratio. The Eddington ratio 
is calculated using the H/3 mass estimate. 
White contours are linearly spaced at in- 
tervals of 50 AGNs per 0.2x0.2-dex bin. 
The BH mass estimates from Mg II are cor- 
related with the Eddington ratio such that 
the mass is underpredicted at low Edding- 
ton ratio and overpredicted at high Ed- 
dington ratio. 
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3.74 
0.13 



3 3.5 4 4.5 



3 3.5 4 4.5 3 3.5 4 4.5 
log FIIHM(H^) 



3 3.5 4 4.5 



Fig. 2.— Distribution of H/3 FWHMs 
for objects having Mg II FWHM within 
the 0.1 dex-wide band shown as the cross- 
hatched bar in each panel. The mean and 
standard deviation in log FWHM(H/?) is 
also indicated for each subsample. Despite 
the overall H/5 and Mg II FWHM distribu- 
tions having similar mean values, the to- 
tal dispersion in log FWHM is 0.18 dex 
for the former and only 0.12 dex for the 
latter, demonstrating that Mg II fails to 
reflect the full range of line widths (and 
hence BH masses). 
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Fig. 3. — Distribution of Eddington ratios 
in bins of bolometric luminosity and red- 
shift, for the redshift range in which Mg II 
estimates of the BH mass contribute to the 
distribution. Dotted histograms employ 
the raw mass estimates from all three emis- 
sion hues. Solid histograms show the result 
after replacing each Mg Il-based mass with 
the corresponding probability distribution 
we derived from H/3. Bolometric luminos- 
ity, redshift, and the number of AGNs are 
shown in each panel. 
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Fig. 4. — Distribution of Eddington ra- 
tios in bins of bolometric luminosity and 
redshift, for objects where both B.f3 and 
Mg II masses were available. Solid his- 
tograms show the distributions for the H/3- 
derived masses, dashed histograms show 
the resultant distributions when our pro- 
cedure is applied to the Mg II masses, 
and the dotted histograms show the raw 
Mg II estimates. The similarity of the 
H/3 (solid) and corrected (dashed) distri- 
butions demonstrates that our method of 
correcting the Mg II masses correctly re- 
covers the true B.f3 Eddington ratio dis- 
tribution. Bolometric luminosity, redshift, 
and the number of AGNs are shown in each 
panel. 
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